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We investigate the physics of coherent quantum phase sUps in two distinct circuits containing 
small Josephson junctions: (i) a single junction embedded in an inductive environment and (ii) a 
long chain of junctions. Starting from the standard Josephson Hamiltonian, the single junction 
circuit can be analyzed using quasi-classical methods; we formulate the conditions under which 
the resulting quasi-charge dynamics is exactly dual to the usual phase dynamics associated with 
Josephson tunneling. For the chain we use the fact that its collective behavior can be characterized 
by one variable: the number m of quantum phase slips present on it. We conclude that the dynamics 
of the conjugate quasi-charge is again exactly dual to the standard phase dynamics of a single 
Josephson junction. In both cases we elucidate the role of the inductance, essential to obtain exact 
duality. These conclusions have profound consequences for the behavior of single junctions and 
chains under microwave irradiation. Since both systems are governed by a model exactly dual 
to the standard resistively and capacitively shunted junction model, we expect the appearance of 
current-Shapiro steps. We numerically calculate the corresponding current- volt age characteristics 
in a wide range of parameters. Our results are of interest in view of a metrological current standard. 

PACS numbers: 74.50. +r, 74.81.Fa, 72.30.+q 



I. INTRODUCTION 

Two physical systems that can be mapped onto each 
other by interchanging the role of position and its canon- 
ically conjugate momentum are said to be related by du- 
ality. If the physical properties of one of the systems are 
known, those of its dual counterpart can be predicted by 
applying the set of duality transformations that accom- 
pany the position and momentum interchange and relate 
the parameters of the two systems. In some special cases, 
duality maps the system onto itself; one then speaks of 
self-duality. An example of a system that exhibits exact 
self-duality is the harmonic oscillator. More frequently 
one encounters systems that exhibit an approximate self- 
duality relating the system's asymptotic behavior in two 
different limiting parameter regimes. 

Duality transformations have been proven useful in a 
variety of situations from a broad range of fields^ includ- 
ing statistical mechanics, condensed matter physics and 
gauge-field theories. Here we will focus on the case of 
Josephson junctions embedded in an electromagnetic en- 
vironment, where the duality associated with the conju- 
gate charge and phase degrees of freedom has be en em- 
ployed to study the circuit's dissipative dynamics'^'^EIIl. 
A Josephson junction, as depicted schematically for a 
current-biased situation in Fig. [l^, is characterized by 
two competing energy scales: the Josephson coupling en- 
ergy Ej = hlc/2e, proportional to the Josephson critical 
current Ic^ and the charging energy Ec = e^/2C, in- 
versely proportional to the junction's capacitance C. A 
self-duality property can be used to relate the junction's 
behavior in the presence of a resistor in the two limit- 
ing cases Ej/Ec > 1 and Ej/Ec < lEEEIIl. This self- 



duality is only approximate, though; as we will discuss in 
more detail below, an exact duality transformation exists 
between the circuit of Fig. [l^ and a different supercon- 
ducting circuit containing a large junction together with 
an additional element: an inductance L^^, see Fig. [l]3. 
This is to be expected somehow: it is customary to de- 
scribe the dynamics of the current-biased circuit depicted 
in Fig. la in terms of a fictitious phase-particle of mass 
C. The dual situation Fig. ^ would then correspond to 
a voltage-biased circuit, the dynamics of which is that of 
a charge-particle of mass L. 

Let us push the duality analysis a little further. The 
phase-particle for the current-biased circuit Fig.[l^ moves 
in a potential U {(/)) which is the sum of a periodic part 
and a linear tilt, see Fig. [2] for an example. The peri- 
odic potential induces the tunneling of Cooper pairs of 
charge 2e, its amplitude Ej = hlc/2e sets the maximum 
supercurrent Ic that can be carried by the junction. The 
tilt of the potential is proportional to the bias current 
lb through the junction. By duality, the charge-particle 
for the voltage-biased circuit Fig. ^ moves in a poten- 
tial U{q) which also has a periodic part and a linear tilt, 
see Fig. [2] Here the periodic part induces a tunneling 
of phase or a "phase slip" whereby the phase difference 
across the junction winds by an amount 27r; accordingly, 
one can speak of a phase-slip junction. A winding of the 
phase with time gives rise to a voltage over the element, 
and the amplitude of the periodic part of the potential 
sets the maximum voltage Vc the phase-slip element can 
sustain. The linear tilt is proportional to the bias voltage 
Vb across the junction. 

Duality thus implies that the I — V characteristics of 
the voltage-biased circuit in Fig.[T]3 can be obtained from 
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FIG. 1: (a) RCSJ-model: Current-biased Josephson junction 
(cross) , bias current h , with Josephson energy Ej shunted by 
a capacitor C and a resistor R. (b) CJRL-model: Voltage- 
biased Josephson junction (cross), bias voltage Vb, shunted by 
a capacitor C and in series with an inductor L and a resistor 
R. 



those of the current-biased circuit of Fig. [T^, by sim- 
ply interchanging the role of current and voltage. This 
has been verified experimentally in Ref.^, where the mea- 
sured I — V characteristics of an underdamped Joseph- 
son junction were found to be exactly dual to those of the 
overdamped Josephson junction measured in Ref.^, in ac- 
cordance with the approximate self-duality exhibited by 
a Josephson junction in these respective limits. 

Probably one of the most important consequences of 
duality in this connection is the case of a junction that 
is irradiated by microwaves (MW) of frequency /. If the 
MW frequency / is commensurate with the frequency 
of the motion of the phase-particle in the periodic po- 
tential, phase- locking occurs yielding Shapiro steps^, in 
the I — V characteristics at well-defined voltages that 
are proportional to multiples n of the applied frequency 
/: Vn = n/i//2e, where h is Planck's constant and e 
the elementary charge. As frequency can be controlled 
with an extreme accuracy, this effect is currently used 
in metrology to define the voltage standard, for a review 
g^ jiQ | ii [ Observation of the dual phenomenon — phase 
locking for the charge-particle yielding Shapiro steps at 
well-defined currents that are multiples of the applied 
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FIG. 2: Example of the tilted washboard potential [/ as a 
function of phase and charge 7vq/e for a phase-particle of 
mass C and a charge-particle of mass L, respectively. The 
dashed line denotes the average tilt, proportional to the bias 
current (phase-particle) through or the bias voltage (charge- 
particle) over the junction. 



frequency 1^ = n2ef — would have far-reaching conse- 
quences for metrology, as this would enable one to define 
a current standard with an unprecedented precision. 

Clearly, the voltage-biased circuit illustrated in Fig.[T]3 
is not the only one dual to that of Fig. [l^. In fact, any 
circuit element that features appropriate phase tunneling 
is a possible candidate for duality. In view of experimen- 
tal implementations and applications it is interesting to 
compare several possibilities. In a recent paper, Mooij 
and Nazaro\^^ proposed exact duality between Fig. [l^ 
and a voltage-biased circuit containing a narrow super- 
conducting wire^^ As we will detail below another 
possibility would be to use a one-dimensional chain of 
Josephson junctions. 

The paper is organized as follows. For pedagogical rea- 
sons and also for the sake of completeness, we start by 
reviewing the physics of the circuits of Fig.[l]in Section[TH 
This enables us to demonstrate the duality principle as 
well as its consequences for superconducting circuits with 

we 



the aid of a relatively simple example. In Section III 



demonstrate that a voltage-biased chain of junctions is 
exactly dual to Fig. [l^. The cur rent- volt age characteris- 
tics of phase-slip junctions both in the absence and in the 
presence of MW irradiation are discussed in Section [IV| 
some perspectives and experimental consequences of our 
theoretical study are discussed in Section [V| 



II. CHARGE-PHASE DUALITY IN CIRCUITS 
CONTAINING A SINGLE JUNCTION 

In this Section we wish to demonstrate that the two su- 
perconducting circuits depicted in Fig. [l] are dual to each 
other. Specifically, we will establish the conditions under 
which this duality holds, treating both circuits entirely 
quantum mechanically. We will discuss the requirements 
to be met such that the classical limit can be taken, lead- 
ing to the usual resistively and capacitively shunted junc- 
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tion (RCSJ) modeF^for Fig.[T^ and its dual counterpart, 
the capacitively shunted junction in series with a resistor 
and an inductance (CJRL) model for Fig. [ij). 



A. Capacitively shunted junction: phase-inertia 

We start our analysis by considering the circuit shown 
in Fig. [l^. It contains a current-biased Josephson junc- 
tion (bias current /fo), shunted by a capacitance C and 
a resistance R. Let (j) be the operator corresponding 
to the phase difference across the junction and Q the 
canonically conjugate charge, such that the commuta- 
tor [Q, ^] = —2ie. The resistor induces dissipation 
that we will account for withi n the framework of the 
Caldeira-Leggett modeW^^^. Hence the circuit pre- 
sented^ in Fig. [l^ can be described by the Hamiltonian 
H — Hq + Hb where 



ffo = ^ + u{4>) , u{4>) 



-Ejcos4>-hIb(j>/2e, (1) 



M. (2) 



i=l 



Here Hb is the Caldeira-Leggett Hamiltonian describ- 
ing a bath of oscillators with frequencies cj^, conjugate 
momenta and positions Pi and X^, the latter coupling 
linearly to the junction's phase operator (j) with coupling 
constants q. 

From the Hamiltonian one can obtain the equation 
of motion for the operator ^, 



nC^/2e + Ic sin (j) = h ^ SI, 



(3) 



where the current SI is related to the momenta of the 
oscillator bath, 



SI . 



(4) 



The bath momenta satisfy the equation of motion 



A+^-A =Q(^. (5) 

Direct integration of Eq. ([5| and substitution of the result 
into Eq. Q yields 



SI . 



j Y{t-t')h4>{t')/2e, 



(6) 



where the first term i{t) is related to the homogeneous 
solution of Eq. ([5|; it is random in nature due to the 
uncertainty with respect to the bath's initial conditions. 
The second term is related to the particular solution of 
Eq. ([5|; it describes the response of the bath to the 



voltage operator ?i(J)/2e through the retarded admittance 
Y(t) with Fourier transform 



f 2e\ cf luo 



\ h J ^ uj (u + iriY — uj ' 

^ ^ i=l f' ^ * 



(7) 



Provided we choose the bath parameters q and coi such 
that 

^Yico)] =(^) nJ2 ^^i^' - ^i) = (8) 

the bath's response is ohmic corresponding to that of a 
resistance R. As a result we can present Eq. ([3| in the 
form 



hC4)/2e + h^/2eR + Ic sin 4) = h 



(9) 



For later use, it is convenient to write Eq. ([9| in a dimen- 
sionless form, dividing both sides by Ic] one then obtains 



(10) 



where r = cOpt with cOp 



{^EjEcY^'^/n the junction's 
plasma frequency, a = {fi/2eIcCY^'^ / R\he dimensionless 

inverse resistance, = Ih/Ic and i = i/Ic- 

Equation (|9| is a non-linear quantum Langevin equa- 
tion, owing its stochastic nature to the presence of the 
random operator i. The statistics of i will be fixed 
by assuming the initial state of the bath to be the 
canonical equilibrium one at temperature T. This, to- 
gether with the harmonic nature of the bath, implies 
that the statistics of the random operator i is Gaus- 
sian with average value (i) =0. In accordance with the 
fluctuation-dissipation theorem, the symmetrized sec- 
ond moment ({z(t), z(0)})/2, where {...,...} denotes the 
anti-commutator, is then characterized by the spectral 
function 



I' 



dte-*({z(t),z(0)})/2 



= nuj^e[Y{u;)] coth{hu;/2kBT) 

= ^coi\i(nuo/2kBT). 
R 



(11) 



A general analysis of the quantum Langevin equation 
([9| is beyond the scope of this article. Here we are in- 
terested in the classical limit of the operator equation 
(|9| where it reduces to the well-known resistively and 
capacitively shunted junction (RCSJ) model, describing 
the classical dynamics of a fictitious phase-particle. The 
capacitor provides the particle's inertia; the correspond- 
ing acceleration is the capacitor's displacement current. 
The resistor provides both velocity-proportional damp- 
ing ~ (j)/R and classical noise i{t). The noise i adds 
to the applied bias current /5 which, together with the 
junction's supercurrent JcSin^ yields the external force 
acting on the particle. 
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A classical interpretation of the operator equation ([9| 
makes sense if we can accurately replace the operators 
by their respective expectation values, (j) ^ (j), i ^ i. 
In particular, we must be allowed to replace (sin^) by 
sin(0) = sin^. For this to be correct, the uncertainty 5(j) 
in the phase must be small compared to the period of the 
sine function. The classical version of ([9| then describes 
the motion of a narrow wave packet of width Scj) <^ 1. 
Such a wave packet can be constructed as a superposi- 
tion of extended phase states, implying an uncertainty 
on the level of the junction charge SQ that exceeds the 
elementary charge e. This means in particular that the 
quasi-classical phase description does not capture effects 
associated with Coulomb blockade^''. 

It is useful to distinguish two cases, according the value 
of the parameter a in Eq. (10): overdamped phase dy- 
namics, corresponding to a > 1 and underdamped phase 
dynamics, a < 1. 

In the overdamped case, the phase dynamics is al- 
ways classical as damping times are naturally short, 
of the order of RC. This is generally achieved in a 
low-resistance environment that avoids charge localiza- 
tion. Classical phase dynamics also requires that the 
environmental noise operator i can be treated classi- 
call}/^^. This corresponds to relatively high tempera- 
tures, such that Siicj) = 2kBT/R, see Eq. ([iT]). Then 
we can replace z by a c-number i such that {%) = and 
(z(t)z(O)) = 2{kBT/R)S{t); the noise is (^-correlated. For 
this to be correct the temperature T should be large 
compared to the characteristic frequency 1/RC. The 
overdamped limit can be analyzed, e.g., by studying the 
Fokker-Planck equation^^ corresponding to the classical 
Langevin equation, as it was done for the overdamped 
case in Ref.!^. 

In the underdamped case, damping times are long and 
wave packet spreading becomes important. This issue 
is particularly relevant when the phase-particle is in the 
running state corresponding to a finite voltage over the 
junction. In the absence of damping, the uncertainty 
S(j) can be kept within limits by applying classical time- 
dependent external forces. The spreading of wave packets 
is governed by a rate proportional to the kinetic energy. 
For the fictitious phase particle this corresponds to the 
charging energy Ec, which is in fact the energy scale 
associated with charge localization leading to Coulomb 
blockade. We require the time-dependence of the exter- 
nal force to be fast on the scale %/Ec- Physically this 
means that we work under conditions where charging ef- 
fects can be ignored. For example, in the case of ex- 
ternally applied microwaves (see Section [lV[) this implies 
their frequency to be larger than Ec/fi- 
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FIG. 3: Lowest two Bloch bands of the spectrum of Hamilto- 
nian Hj^ taking the ratio Ej/Ec — 2. Note that the bands 
are separated by an energy of the order of the plasma fre- 
quency, Tiujp — 4:Ec' 



inductive and resistive series environment (inductance L 
and resistance R^^. Its Hamiltonian is given by H' = 
Hj + Hl with 



Here H^^ 



Hj = 



91 
2C 



Ej cos ; Hl = 



(06 -^R- 0? 

2L 



(12) 

where 6 is the phase difference across the junction, con- 
jugate to the charge Q, such that [Q^O] = —2ie. The 
combination 06 — ^i? — ^ is the phase difference across the 
inductance. It contains the total phase difference across 
the circuit 05, which is an external parameter determined 
by the applied voltage bias 14, such that 05 = 2eVb/h. 
The operator (j)R is the phase difference across the resis- 
tor. We express it in terms of the oscillator bath positions 
according to the relation 



(13) 



with coupling constants A^. The dynamics of the bath 
degrees of freedom, accounting for dissipation due to the 
resistor, is governed by the Hamiltonian H^^ 



2 



(14) 



B. Junction in series with an inductance: 
charge-inertia 

Next consider the circuit depicted in Fig. [T]3, where 
a single Josephson junction is embedded in a combined 



where 11^ are the bath momenta conjugate to S^; Coi are 
the bath frequencies. 

Let us first consider the junction Hamiltonian Hj^ 
which describes a particle of mass C moving in a periodic 
potential. The spectrum of this Hamiltonian consists of 
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Bloch bandPI^, see Fig. |3| W e can use the so-called 
quasi-charge representatiorP^^ and write Hj = eo(^), 
assuming the junction's dynamics to be restricted to 
the lowest Bloch band eo{q). The commutation rela- 
tion of the quasi-charge operator q and the phase 6 reads 
[q^O] = —2ie. The physical meaning of the quasi-charge 
becomes clear if we consider its dynamics, which is gov- 
erned by the equation of motion 



2eL^ 



2e dH' 
ft dSh ' 



(15) 



Since 0^ is the externally fixed phase drop over the entire 
circuit, the last term on the right-hand-side corresponds 
by definition to the total current fiowing through the cir- 
cuit. We next calculate = i[H' ^0]/fi and find 



2e dcQ^q) 



(16) 



provided that interband transitions can be ignorecP^. 
Then the operator deo{q)/dq corresponds to the voltage 
drop over the junction. A straightforward calculation of 
q now yields the equation of motion for the quasi-charge. 
We find 



Lq^deo{q)/dq = Vb^6V, 
where the operator 

h 



SV 



2e 



(17) 



(18) 



The bath positions satisfy the equation of motion 

Ei^Co^Ei = Xih4/2e. (19) 



Direct integration of Eq. (19) yields the sum of the ho- 



mogeneous and the particular solution, substitution into 
Eq. (pl yields 



5V = v{t)- J Z{t-t')'q{t'). 



(20) 



Similar to the situation discussed above for the operator 
5/, the first term v{t) is the voltage noise related to the 
uncertainty on the level of the initial conditions for the 
homogeneous solution of Eq. ( 19 ). The second term is re- 



lated to the particular solution of Eq. ( 19 ) and describes 



the response of the bath to an applied current q through 
the retarded impedance Z{t). The Fourier transform of 
the latter is 



Z{lo) 



2 oo 



~,2 • 



(21) 



If we choose the bath parameters and Coi such that 



2 OO 



Coi) = R. (22) 



the bath's response is ohmic corresponding to that of a 
resistance R. Assuming the bath's initial state to be a 
canonical equilibrium one, we find that the voltage noise 
is characterized by the two lowest cumulants, {v) = 
and ({'O(t), '0(0)})/2, the latter satisfying the fiuctuation- 
dissipation theorem. 



S.ico) = j dte'^'{{v{t),v{0)})/'^ 

= hu^e[Z{ij)] coth{hu/2kBT) 

= hujRcoth{huj/2kBT). (23) 

Now consider the limit Ej Eq^ where eo{q) = 
— Ao cosnq/e corresponds to a purely sinusoidal band in 
quasi-charge representation with a bandwidth given 



Ao = 16^/EjE^{Ej/2Ec)'^^e-V^^^. (24) 



Then, together with the choice ([22| for the oscillator bath 
parameters, Eq. (17) takes the form 



Lq-\- Rq-\- Vc sin nq/e = Vi} 



(25) 



where Vc = 7rAo/e is the critical voltage. If both sides of 
Eq. (25) are divided by Vc^ we obtain the dimensionless 
form 



d^q/dr'^ + pdq/dr + sing = Vb 



- V. 



(26) 



Here, r' = coct with coc = {TrVc/eLy^'^ dual to the plasma 
frequency, q = nq/e^ the dimensionless resistance p = 
R{e/7rVcL)^/^ H = VhlVc and S = vlV^. 



Comparing the quantum Langevin equation Eq. (25) 



with the corresponding one for the phase dynamics, 
Eq. ([9|, we conclude that they are exactly dual to each 
other. In other words, Eqs. (25) and (l9| map onto each 



other when exchanging the role of quasi-charge q and 
phase (/), such that nq/e (/>, accompanied by the dual- 
ity transformations e/n ^ h/2e, L ^ R ^ 1/R, and 
V ^ I. This establishes the exact duality between the 
two circuits. Fig. [l^ and b. 

Equation ( [25| ) has a simple physical interpretation in 
the classical limit: it describes a fictitious charge-particle 
with inertia L, provided by the inductor. The sum 
of the applied bias voltage V5 and the resistor-induced 
noise v drops over the series configuration formed by 
the junction, the inductor and the resistor. The charge- 
particle moves in a tilted washboard potential U{q) = 
—So cosnq/e — V^q and experiences velocity-proportional 
damping Rq. 

As for its dual counterpart, the classical interpretation 
of the operator equation ( 25 ) hinges on the replacement 
of the operators q and v by their respective expectation 
values q and v. This means in particular that we must 
be allowed to replace {smnq/e) by sin(7rg/e) = smirq/e. 
For this to be correct, the uncertainty in the charge Sq 
must be small compared to the elementary charge e. The 
classical version of (25) then describes the motion of a 



narrow wave packet of width 5q <^ e. In phase rep- 
resentation this is consistent with the reahzation of an 
extended Bloch state. 

As for the case of phase dynamics, it is useful to distin- 
guish two situations, according the value of the param- 
eter p in Eq. (26): overdamped charge dynamics, corre- 



sponding to p > 1 and underdamped charge dynamics, 

p<i. 

In the overdamped case, the charge dynamics is always 
classical as damping times are naturally short, of the or- 
der of L/R. Note that L here constitutes an additional 
element in addition to the junction, unlike the case of 
phase dynamics discussed above where the capacitance 
C is a property intrinsic to the junction. Overdamped 
charge dynamics is generally achieved in a high-resistance 
environment, which favors charge localization. Classi- 
cal charge dynamics also requires that the environmental 
noise operator v can be treated classicall}'^. This implies 
working in the high-temperature limit ksT > hR/L^ 
such that V can be replaced by a c-number with (v) = 
and {v{t)v{0)) = 2kBTR5(t) in accordance with Eq. (|23|. 
The resulting classical Langevin equation can be an- 
aly zed t hrough the corresponding Fokker-Planck equa- 
tior]IHE3. 

In the underdamped case, damping times are long and 
wave packet spreading becomes important. This issue is 
particularly important when the charge-particle is in the 
running state corresponding to a finite current through 
the junction. In the absence of damping, the uncertainty 
5q can be kept within limits by applying classical time- 
dependent external forces. The spreading of wave packets 
is governed by a rate proportional to the kinetic energy. 
For the fictitious charge particle this corresponds to the 
inductive energy El = ^o/2^. which is in fact the en- 
ergy scale relevant for phase localization related to the 
Josephson effect. We require the time-dependence of the 
external force to be fast on the scale ?i/El- Physically 
this means that we work under conditions where phase 
slip events are not suppressed. For example, in the case 
of microwaves discussed in Section [IV| below, this implies 
their frequency to be larger than El/TI' 

We conclude this section by summarizing the condi- 
tions under which the above exact duality is obtained. 
First of all, we assume the quasi-charge dynamics of 
the voltage-biased circuit Fig. [TJd to be determined by 
the lowest Bloch band only. The lowest band is sepa- 
rated from the next one by the plasma frequency fiuOp = 
\/^EjEc' Since we do not consider here inter-band tran- 
sitions, we thus assume all energies to be smaller than a;^. 
Second, in order for the Bloch band to be purely sinu- 
soidal, we need to impose the condition Ej ^ Eq for the 
voltage-biased junction of Fig. [TJd. 



III. JOSEPHSON JUNCTION CHAIN 

In the previous section we have found that a voltage- 
biased Josephson junction with a large ratio Ej/Ec in se- 
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FIG. 4: Josephson junction chain. 



ries with an inductance constitutes a phase-slip element. 
We also saw that the observation of well-defined quasi- 
charge dynamics requires charge-fiuctuations very much 
smaller than a single charge e, a condition that can be 
obtained by a large inductance near the Josephson junc- 
tion. Experimentally it is not so simple to realize large 
magnetic inductances very close to the sample. An al- 
ternative is to fabricate a large kinetic inductance using 
a superconductor, and in particular a Josephson junc- 
tion chain^^^^^. It seems therefore quite natural to an- 
alyze the possibility to realize a phase-slip junction from 
a Josephson junction chain, which we investigate in this 
section. The central idea is that the phase-slip itself oc- 
curs on only one of the junctions of the chain; the phases 
on the other junctions just perform small Josephson os- 
cillations, thereby providing the necessary inductance. In 
order to demonstrate this idea, we will closely follow the 
paper by Matveev et al.^, who studied the low-energy 
properties of Josephson junction chains in terms of quan- 
tum phase slips. 

We start our analysis by considering the Josephson 
junction chain depicted in Fig. [4j It is a series arrange- 
ment of N identical junctions, each with Josephson en- 
ergy Ej = hlc/2e and charging energy Ec = e^/2C Let 
Qi be the charge on the i-th junction and Oi the conju- 
gate phase difference. We consider the nearest-neighbor- 
capacitance limit and assume the chain to be subjected 
to an external phase (f). The Hamiltonian can then be 
written as 



N 



N 



1=1 i=l 

(27) 

Ignoring the charging energy for the moment, we find 
the classical ground state configuration of the chain's 
phases Oi by minimizing the Josephson coupling energy, 
thereby satisfying the constraint. The corresponding 
configuration can be easily found in the limit N ^ 1 
and is given by Oi = (j)/N. The resulting Josephson en- 
ergy hence reads = Ejcj)^ /2N. This is the inductive 
energy stored in the chain; the corresponding effective 
inductance is given by Lar = TiN/2eIc' 

Now consider a phase slip event occurring on one of the 
junctions, say the jth junction, such that Oj Oj + 27r. 
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Due to its periodicity as a function of the Joseph- 
son energy of the junction j does not change. How- 
ever, the constraint = ^ is violated after such a 
phase-shp event. In order to accommodate the phase- 
shp without violating the constraint, the phase differ- 
ences 6i over the other junctions change slightly, from 
(j)/N to {(j) — 27t)/N. Correspondingly, the Josephson 
energy of the chain changes from Eq = Ejcj)^ /2N to 
El = Ej(y(j) — 2itY /2N . In a similar way one shows that 
the classical energy needed to accommodate m phase- 
slips without violating the constraint is given by Ej^ = 
Ej{(j) — 27rm)^/2A/'. We thus conclude that the ground 
state of the chain generally is one that contains a fixed 
number of phase slips for almost any value of the exter- 
nal phase ^, except for the special values (j) = 7r(2m + 1) 
where the energies Em and Em-\-i are degenerate. Quan- 
tum fiuctuations induced by the small but finite charging 
energy Ec lift this degeneracy: they give rise to a non- 
vanishing amplitude Aq for a phase-slip event to occur. 
We denote the state of the chain with m phase slips by 
\m). Taking into account the fact that a phase slip can 
take place on any of the N junctions, we can write the 
total Hamiltonian for the chain as 



^^[|m+l)(m| + /i.c.]. (28) 



Next introduce the operator conjugate to the phase- 
slip number m. Specifically, m] = — ze/yr such that 
the operator e^^^l^ is a raising operator with e^^^l^\m) = 
|m + 1). Using this representation, the Hamiltonian be- 
comes 



^ar = (^j/27V)(2^m- 



NAQCO^irq/e. (29) 



A physical interpretation of the operator q can be 
obtained by calculating ^; the result reads q = 
{2e/h)dH /d(p^ which by definition is the operator corre- 
sponding to the total current through the chain. Hence, 
in analogy with the result (15), q can be interpreted as 
the global charge of the chain. It is easy to see that 



Hamiltonian Eq. (29) has the same form as the quasi- 



charge representation of Hamiltonian Hq for a single 
junction in series with an inductor, presented in the pre- 
vious section. It is therefore straightforward to analyze 
the case of a voltage-biased chain embedded in a resis- 
tive series environment, repeating the steps presented in 
Section [HB) One finds that the quasi-charge dynamics is 
governed by the equation 



Lar<7 Rq^ K,ar siu TTq/c = 



where 

4>b = 
sented 

^c,ar 
P = 



2eVb/h. 
in the 

R{e/nVc 



NVc and we used the fact 
This result can also be 
dimensionless form ( 26 ) , 



(30) 

that 
pre- 

with 
and 



Cariar)^/' = (i?/iV) (26^ 1,^^) ■ 

This result refiects the intuitive argument discussed at 
the beginning of this section: the inertia Lar of the charge 



dynamics is provided by the chain itself. Moreover, it is 
given by N times the nonlinear inductance h/2elc of a 
single junction in the chain. This means that Lar can 
be tuned in principle, either by tuning N or by using 
SQUID loops instead of single junctions as in'^, such 
that Ic can be tuned with a magnetic fiux. We also note 
that the critical voltage of a Josephson junction chain 
is N times larger than the one of a single Josephson 
junction. This is relevant for the discussion in the next 
section, where we will analyze the I-V characteristic of 
a phase-slip junction under microwave irradiation. As 
the width of the appearing current steps scales with the 
critical voltage of the phase-slip junction, a Josephson 
junction chain has necessarily larger current-steps that 
are as a consequence more robust against voltage noise. 



IV. PHASE-SLIP JUNCTION UNDER THE 
INFLUENCE OF MICROWAVE IRRADIATION 

As it was already mentioned in the Introduction, it 
is of interest to study the behavior of phase-slip junc- 
tions under the infiuence of microwaves. In view of du- 
ality, we expect the current- voltage characteristics to ex- 
hibit steps at well-defined values of the current that are 
multiples of the microwave frequency, so-called current 
Shapiro steps. In order to demonstrate this, we have nu- 
merically integrated Eq. ( [25| ) in the classical limit. We 
set Vb{t) = V -\- Fmw sin(c3MW^), and ignore the effect 
of fiuctuations. The results are plotted in Figs. [5] and [6] 
for two choices of the dimensionless damping parameter 
p, corresponding to the overdamped and underdamped 
limit, respectively. 

Figures[5^ and|6^ show the current- voltage characteris- 
tics in the absence of microwaves, Vmw = 0. In the over- 
damped case, Fig.|5^, the current remains zero as long as 
the voltage V is smaller than the critical voltage Vc'. the 
phase-slip junction is in the Coulomb blockade regime. 
Beyond Vc^ the current rapidly increases until it reaches 
the linear asymptote I ^ V/R^ where the junction is in 
the superconducting state and the voltage drops entirely 
over the resistor. This characteristic corresponds in fact 
to the so-called Bloch nose, which has been studied in the 
experimeniP in the presence of thermal fiuctuations. It 
was found that such fiuctuations induce a smooth inter- 
polation between the Coulomb blockade and the super- 
conducting branch, in agreement with theor}!^. In the 
underdamped situation. Fig. |6^, the two branches in the 
characteristic co-exist in the Coulomb blockade region. 
Since both branches are accessible, we generally expect 
to find hysteretic behavior of the phase-slip junction. To 
the best of our knowledge, this limit has not yet been 
accessed experimentally. The behavior shown in Figs. [5^ 
and[6^ have a dual counterpart in the usual RCSJ-model, 
where both the overdamped and the underdamped limit 
have been studied in experimentJ^l^. 

If microwaves are applied, steps appear at well-defined 
values of the current that are multiples of the applied 
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FIG. 5: Current-voltage characteristics for an overdamped 
phase-slip junction with p = 2, (a) without MW irradiation, 
(b) under MW-irradiation with amplitude Vmw = 5^ and 
frequency umw — l.ScJc- 
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FIG. 6: Cur rent- volt age characteristics for an underdamperd 
phase-slip junction with p = 0.15, (a) without MW irradia- 
tion, (b) under MW-irradiation with amplitude Vmw = 314 
and frequency cjmw = Scjc- 



V. CONCLUSIONS AND OUTLOOK 



MW frequency, as can be seen in Figs, ^jp and ^jp. The 
steps are clearly visible; they are integrated within the 
overall I-V characteristic, the average slope of which 
remains determined by the resistance R. In the over- 
damped case, there is a one-to-one correspondence be- 
tween current and voltage, as for the case without MW. 
In fact, the shape of each step in Fig. [sJd appears to 
be a replica of the characteristic in the absence of MW, 
Fig. [5^. The one-to-one correspondence between volt- 
age and current found in the overdamped limit is lost 
in the underdamped case: various current steps appear 
within the same voltage interval. In order to obtain the 
result shown in Fig. |6fc), Eq. (25) had to be integrated 



for a range of initial conditions on the quasi-charge for 
each value of the DC voltage V. Note that the resistive 
branch is absent in the regions of overlapping steps; this 
is of interest as it possibly makes it easier to lock on a 
given step in the experiment. 



We have studied two circuits, a single Josephson junc- 
tion in an inductive-resistive environment and a Joseph- 
son junction chain, in view of the realization of the dual 
of the Josephson effect. In both cases we elucidate the 
importance of the role of the inductance in order to re- 
duce charge fluctuations. Duality between these two cir- 
cuits and a single Josephson junction is valid in the quan- 
tum and classical regime. Here, we derived I{V) charac- 
teristics in the classical quasi-charge regime. In case of 
larger quasicharge fluctuations going beyong the classi- 
cal regime, perturbation theory can be applied in analogy 
to the P{E) th eory in Josephson junctions governed by 
phase-dynamic j^ '^^'^^i 

Until now there are only a few e xperiments dealing 
with quasi-charge dynamic j^ l ^^ * ^^ * ^^ . We believe that 
both circuits are experimentally feasible and they are 
of particular interest in terms of the realization of cur- 
rent Shapiro steps. The successful realization of such 
an experiment would link the frequency to the current 
by a quantum electrical recipe and close the metrologi- 
cal triangle. Ultimately, the quantum metrological trian- 
gle experiment would enable a consistency check of the 
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fundamental constants of nature, the electron charge e 
and Planck's constant /i, and as a consequence link the 
kilogram (the only unit still defined on the basis of one 
artefact, prototype of the mass kept in metrology insti- 
tutions), to the Planck constant. The circuit based on 
a single Josephson chain seems from our point of view 
the most promising as it enables to realize large current 
steps. 
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